Astronomy & Astrophysics manuscript no. HD168112 


February 2, 2008 


(DOI: will be inserted by hand later) 





Non-thermal radio emission from O-type stars 



I. HD 168112 



in 
o 
o 



R. Blomme 1 , S. Van Loo 1 , M. De Becker 2 , G. Rauw 2 *, M.C. Runacres 1 , D. Y. A. Setia Gunawan 3 , and J.M. 

Chapman 3 

1 Royal Observatory of Belgium, Ringlaan 3, B-l 180 Brussel, Belgium 

2 Institut d' Astrophysique, Universite de Liege, Allee du 6 Aout, Bat B5c, B-4000 Liege (Sart-Tilman), Belgium 

3 Australia Telescope National Facility, PO Box 76, Epping, NSW 2121, Australia 

Received date; accepted date 



(N 

> 
in 

(N 
O 

in 
o 

Of 

6 



Abstract. We present a radio lightcurve of the 05.5 111(F) star HD 1681 12, based on archive data from the Very Large Array 
(VLA) and the Australia Telescope Compact Array (ATCA). The fluxes show considerable variability and a negative spectral 
index, thereby confirming that HD 168112 is a non-thermal radio emitter. The non-thermal radio emission is believed to be 
due to synchrotron radiation from relativistic electrons that have been Fermi accelerated in shocks. For HD 168112, it is not 
known whether these shocks are due to a wind-wind collision in a binary system or to the intrinsic instability of the stellar 
wind driving mechanism. Assuming HD 168112 to be a single star, our synchrotron model shows that the velocity jump of the 
shocks should be very high, or there should be a very large number of shocks in the wind. Neither of these is compatible with 
time-dependent hydrodynamical calculations of O star winds. If, on the other hand, we assume that HD 168112 is a binary, 
the high velocity jump is easily explained by ascribing it to the wind-wind collision. By further assuming the star to be an 
eccentric binary, we can explain the observed radio variability by the colliding-wind region moving in and out of the region 
where free-free absorption is important. The radio data presented here show that the binary has a period of between one and 
two years. By combining the radio data with X-ray data, we find that the most likely period is ~ 1.4 yr. 
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1. Introduction 

Radio emission from most hot stars is due to the thermal free- 
free emission by material in their stellar winds. A significant 
fraction of these stars, however, also show evidence of non- 
thermal emission. The radio fluxes of non-thermal sources are 
characterized by a zero or negative spectral index 1 , signifi- 
cant variability and a high brightness temperature (Bieging 
et al. 119891 . White J1985i attributed the non-thermal radia- 
tion to relativistic electrons that spiral in a magnetic field and 
thereby emit synchrotron radiation. These electrons are accel- 
erated by the first-order Fermi mechanism (Bell 1978 1 in the 
shocks present in the stellar wind. In a binary, the relativistic 
electrons are accelerated in both shocks formed by the collid- 
ing winds (Eichler & Usov 1993 1. In a single star, the shocks 
are due to the instability of the radiative driving mechanism 
(e.g. Owocki et al.[Tp88 1. 
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For Wolf-Rayet stars, Dougherty & Williams (2000 1 
showed that non-thermal emission is strongly correlated with 
binarity: 7 out of the 9 non-thermal WR emitters are identi- 
fied as binaries. Not all WR binaries are non-thermal emitters: 
if the period is sufficiently short (P 5; 1 yr), the synchrotron 
emission is completely absorbed by the high free-free opacity 
in these stellar winds. 

For O stars, the link between non-thermal emission and 
binarity is less clear. While colliding winds are an explana- 
tion in some cases, there remain a number of apparently sin- 
gle O stars that are non-thermal emitters. Relativistic electrons 
in these stars would then be accelerated in wind-embedded 
shocks. Time-dependent hydrodynamical calculations confirm 
that the instability in the radiative driving mechanism results in 
a significant number of shocks (e.g. Runacres & Owocki 2005 
and references therein). Models for non-thermal emission of O- 
type stars, based on this mechanism, have been developed by 
Chen & White JT9941 and Van Loo et al. J20041 12005b . 

The purpose of the present paper is to study the non-thermal 
radio emission of the apparently single star HD 1681 12. We do 
this by collecting all available archive data. We explore both the 
single-star and binary scenario, to see which one best explains 
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Table 1. The parameters of HD 1681 12 from the SIMBAD cat- 
alogue, Leitherer dl988> . Bieging et al. (1989 1 and Conti & 
Ebbets (tT9T7l . 
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the characteristics of the non-thermal emission. The main char- 
acteristics we study are the radio flux levels, spectral shape and 
time-scale of the variability. 

HD 1681 12 belongs to the Ser OB 2 association (Cappa et 
al. l2002> and lies inside NGC 6604, located at the core of the 
H II region S54 (Georgelin et al. 1973 1. The stellar parameters 
and wind parameters used in this paper are listed in Tabled De 
Becker et al. ( 2004 1 classified HD 1681 12 as 05.5 III(f + ). They 
also showed that optical spectra do not reveal radial velocity 
variations: thus, HD 168112 is not a binary, or it has a low 
inclination orbit or a long period. 

Abbott et al. ( 1985} were the first to identify HD 168112 
as a non-thermal emitter. Radio observations indicate a zero 
or negative spectral index as well as considerable variability. 
Moreover, even at minimum, its 6 cm radio flux is a factor 
of 10 higher than that expected from free-free radiation of the 
mass loss rate as determined from Ha, which is a further indi- 
cator of non-thermal emission (Bieging et al. 1989 De Becker 
et al. l2004l . 

De Becker et al. presented two radio observations sepa- 
rated by 6 months, showing a factor of 5-7 increase. Quasi- 
simultaneous XMM-Newton observations show anti-correlated 
behaviour: there is a decrease of ~ 30 % in the X-ray flux. De 
Becker et al. interpret such behaviour as pointing to a binary 
scenario, where a significant fraction of the X-ray emission is 
produced by the hot plasma heated by the wind-wind interac- 
tion. The present paper discusses a much larger number of radio 
observations. These will allow us to quantify the non-thermal 
properties of HD 1681 12 and thus to better evaluate the merits 
of the single-star versus binary scenario. 

The remainder of this paper is structured as follows. In 
Sect. |2 we present the data. In Sect. we explore both the 
single-star and binary scenario and in Sect. |4] we present our 
conclusions. 



2. Data 

Tables IA.1I and IA.2I list the radio observations centred on, 
or near to, HD 168112. These data were collected from the 



archives of the NRAO Very Large Array 2 (VLA) and the 
Australia Telescope Compact Array 3 (ATCA). Many of these 
observations have not been published previously. To avoid 
introducing systematic shifts between different data sets, we 
decided to re-reduce all observations in a consistent way. 
Details of the reduction and the measured fluxes are given in 
AppendixlAl 

2.1. Long-term variability 

The VLA and ATCA fluxes at 2, 3.6, 6 and 20 cm are plotted 
as a function of time in Fig.^ Not all 20 cm upper limits were 
included in this figure, only those which are lower than 13 mJy. 
The single 13 cm detection is not plotted. The 90 cm observa- 
tions are not shown as they only provide upper limits. Note that 
there may be slight differences in wavelength for observations 
we classed in the same wavelength band. The exact frequencies 
for each observation are listed in Table lA.31 

The detections at 2, 3.6 and 6 cm in Fig^dearly show the 
variability that is an identifying characteristic of non-thermal 
emission. Observations differing by more than 0.4 yr show sig- 
nificantly different fluxes. At certain times, HD 168112 is in a 
low state (around year 1992.4 and 2002.2), while at other times 
it is in a high state (1985.0, 1993. 1 and 2001 .8). The flux differs 
by a factor of ~ 20 between these states. At 1984.2 and 2002.9, 
the fluxes are intermediate between the two states. Although 
few detections have been made at 20 cm, there are some defin- 
ing a high state (1985.0, 2001.8) and some of the upper lim- 
its are sufficiently sharp to confirm the presence of a low state 
(2002.2). 

The flux variations are not due to the fact that we are resolv- 
ing the wind of the star in the highest-resolution configurations. 
A high-resolution observation typically has a beam of 075 di- 
ameter, corresponding to ~ 13 000 R* at 2 kpc distance. This 
value is much higher than the expected size of a synchrotron 
emitting region, either in the single-star or binary scenario. We 
also note that the 20 cm images contain another source (GPS 
G018.565+01.756,at ~ 11' from HD 168112): the flux of this 
strong source (~ 50 mJy) varies by only 20 % between the 
HD 1681 12 high (1985.0) and low (2002.2) state. We therefore 
consider that an instrumental effect is unlikely to be the expla- 
nation for the observed variability. 

FigureQJalso shows that the spectral index between 3.6 (or 
2) and 6 cm is negative in most cases. This is a further indica- 
tion that the emission is non-thermal. The only exception is at 
2002.2, where there is a large error bar that includes the ther- 
mal value {a = +0.6). This does not mean, however, that we 
are observing the underlying free-free emission of the wind. 
As De Becker et al. ( 2004 1 already pointed out, the expected 
thermal flux is only 0.03 mJy at 6 cm, which is a factor of 10 
lower than the observed flux. We further note that the most ex- 



2 The National Radio Astronomy Observatory is a facility of the 
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eration as a National Facility managed by CSIRO. 
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Fig. 1. The 2, 3.6, 6 and 20 cm fluxes of the VLA and ATCA archive data, plotted as a function of time. The bottom panel shows 
the spectral index between 3.6 and 6 cm (or 2 and 6 cm if there is no 3.6 cm observation). Observations are labelled by their 
programme name. Only those 20 cm upper limits below 13 mJy have been included. 



treme values of the spectral index occur when the most extreme 
values of the flux are reached. When the flux is high, we find 
a very negative spectral index (a - -1.1 + 0.4 at 2001.1 and 
a = -0.78+0.08 at 1985.0). When the fluxes are low (2002.2), 
we find a = +0.6 + 0.6. 

The spectral index between 6 and 20 cm provides little ad- 
ditional information and is therefore not shown in Fig. [2 It is 
flatter than the 2-6 cm one at 1985.0. Around 2002.9, the 6- 
20 cm index is positive, indicating a turnover in the spectrum. 



Other observations do not allow a significant determination of 
the 6-20 cm spectral index. 

2.2. Short-term variability 

The shortest time-scales are covered by the observations with 
programme name AB671 (around 1993.1). These data consist 
of 6 observations at both 3.6 and 6 cm, spread roughly uni- 
formly over ~ 0.1 yr (1993 January 21 to February 14). Figure^] 
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Fig. 2. The 3.6 and 6 cm fluxes of the AB67 1 observations. The 
solid line is the flux measured on the combined dataset. 

shows that observations No. 2 and 4 deviate from the average 
by 5-10 %. While this may suggest variability, it will be shown 
that this variability is spurious. 

A typical test to detect an instrumental cause for appar- 
ent variability is to measure another source on the same im- 
age and see whether any flux changes are correlated to the ob- 
ject changes. Unfortunately, only one additional source is vis- 
ible on the 3.6 and 6 cm images (at RA(J2000)=18 h 18 m 40 s , 
Dec(J2000)=-12°07') and it is too weak (~ 0.5 mJy) to see if 
its flux changes are correlated with those of HD 1681 12. 

A possible cause of apparent variability is the presence 
of large atmospheric phase changes that are not sufficiently 
well corrected for. In the reduction of the AB67 1 dataset, self- 
calibration was applied to follow the phase changes on a shorter 
time-scale than the phase calibrators allow (Appendix IA.2> . 
Because observation No. 4 had the lowest flux in the stan- 
dard reduction (i.e. without self-calibration), we suspect that in 
this specific case significant phase changes occurred on shorter 
time-scales than could be corrected for with self-calibration. 
The fact that the flux for No. 4 is low is therefore not signifi- 
cant. 

Observation No. 2 deviates from the solid line in Fig. |2] 
by 1.33 <t at 3.6 cm and 1.64 <x at 6 cm. A value of 1.64 <x 
corresponds to a 10 % (two-sided) probability of a Gaussian 
distribution. As we have 6 observations, it is not significant that 
one of them deviates by 1.64 cr. We therefore conclude that 
we find no evidence for significant variability in the AB671 
dataset. This allows us to combine all AB671 data into a single 
dataset and apply our reduction technique to it. The result for 



both wavelengths is given in Table |A~T1 and shown as a solid 
line in Fig. [2] 

Some of the other data allow us to look for variability on 
the 0.1-0.2 yr time-scale. The two 6 cm observations AA28 
and AA29 (at 1984.2) are ~ 0.1 yr apart and the flux is con- 
stant, within the error bars. There are three AC1 16 20 cm ob- 
servations (around 1985), which cover a ~ 0.2 yr time-scale. 
Unfortunately, one of the observations only gives a (high) up- 
per limit and another has a large error bar. No significant infor- 
mation about variability can be derived from these data. 

We do find clear variability on the time-scale of 0.45 yr 
between the C978 (2001.8) and the AB1048 (2002.2) observa- 
tion. The flux drops from the high state to the low state at 3.6, 6 
and 20 cm. This provides the shortest time-scale on which the 
radio observations show variability. 

3. Interpretation 

3. 1 . Single-star scenario 

The non-thermal radiation is due to synchrotron emission by 
relativistic electrons that are Fermi accelerated in shocks. In 
the single-star scenario, the shocks are wind-embedded (White 
1985 1 as a consequence of the instability of the radiative driv- 
ing mechanism. Time-dependent hydrodynamical calculations 
confirm that the wind is pervaded by shocks that can survive 
to large distances (e.g. Runacres & Owocki 2005, and ref- 
erences therein). Models for non-thermal emission of O-type 
stars, based on the White (1985 1 mechanism, have been de- 
veloped by Chen & White (fT994l l and Van Loo et al. (2004 
2005 1. These models are indeed able to explain the observed 
non-thermal radio fluxes of Cyg OB2 No. 9 and 9 Sgr, though 
Van Loo et al. |2005t raise some questions about the high shock 
velocity jump required for their fit of the Cyg OB 2 No. 9 data. 

In this section, we apply the Van Loo et al. (2005 1 model 
to HD 168112. While shocks have a large range of compres- 
sion ratios and velocity jumps, Van Loo et al. showed that the 
strongest shocks dominate the synchrotron emission. For sim- 
plicity, the model assumes a single shock to be present in the 
wind. (This assumption is of minor importance as the results 
can be re-interpreted in terms of multiple shocks - see below.) 
The model is then described by only four parameters: the posi- 
tion of the shock R$ , the compression ratio x an d the velocity 
jump Am of the shock and the surface magnetic field B„. 

In the model, the momentum distribution of the relativis- 
tic electrons is calculated taking into account the electron ac- 
celeration by the Fermi mechanism (Bell I1978> and the de- 
celeration due to inverse-Compton and adiabatic cooling. The 
synchrotron-emitting relativistic electrons form a narrow layer 
downstream of the shock. The synchrotron flux is then calcu- 
lated (including the Razin effect) and corrected for free-free ab- 
sorption in the wind, giving the predicted radio flux. We assume 
that for this star, v m = 250 km s _I (similar to White (1983 f° r 
Cyg OB 2 No. 9) and that the wind temperature is 0.3 times the 
effective temperature (Drew 1989 1. For further details of the 
model, see Van Loo et al. J2005> . 

By calculating a grid of models, we can see which combi- 
nations of parameters explain the observations. To constrain the 
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Fig. 3. (a) Left panel: The combinations of %, Rs and B* that fit the AO 16 (1984 December 21) observation for HD 168112 lie 
in the grey volume on the figure. Projections on two planes are plotted to situate the solutions in the parameter space, (b) Right 
panel: The values for the surface magnetic field B* and the shock velocity jump Am are shown for all possible combinations of 
the model parameters that explain the radio observations of HD 1681 12. 



fit parameters sufficiently well, it is necessary to have (nearly) 
simultaneous radio observations at three or more wavelengths. 
We applied the model to the AC 116 (1984 December 21) ob- 
servation, which has detections at three wavelengths (2, 6 and 
20 cm). Note that this observation is representative of the high 
flux state of HD 168112. We calculated the fluxes at the three 
wavelengths for a grid of Rs, x> Am ar, d B*. The star and 
wind parameters needed in the model are given in Tabled All 
combinations of the model parameters that fit the observations 
within the error bars are shown in Fig. [3] 

Figure |5Ji shows that only moderately strong to strong 
shocks (x - 2.7 - 4), located between 220 R* and 500 R* 
can explain the observations. The surface magnetic field is 
less well constrained: we find results for B, from 30 G up- 
ward, to values which are above the generally accepted de- 
tection limit for O stars (of order 100 G, Mathys 1999). The 
most important result, however, is that the shock velocity jump 
is 600 - 1800 kms -1 (Fig. OJ?). Such values are considerably 
higher than those found from time-dependent hydrodynamical 
calculations: at 200 - 500 R„, Au is typically only 50 kms -1 
(Runacres & Owocki l2005l . 

One can lower the required Au by re-interpreting the results 
of this single shock model in terms of multiple shocks. Van Loo 
et al. (2005 1 showed that multiple shocks also fit the observa- 
tions provided that the sum of their Au 3 equals the single shock 
Am 3 (all other parameters being the same). To attain the reason- 
able value of Au = 50 kms 1 would require more than 1 700 
shocks. We cannot put these shocks at arbitrary positions in the 
wind (Van Loo 2005 1. In the single-shock model, the resulting 
20 cm flux is determined by both the intrinsic synchrotron ra- 
diation and the free-free absorption. If we put the shocks in a 
region well beyond the location we found for the single shock, 
the free-free absorption will be less effective and the 20 cm flux 
will become higher than observed. We therefore have to put the 
1 700 shocks in a region near the location we found for the sin- 



gle shock. From hydrodynamical calculations, it is obvious that 
we cannot put that many shocks in so small a region. 

Another observation suitable for the Van Loo et al. ( 2005 1 
model is the C978 one, which provides simultaneous observa- 
tions at 3.6, 6, 13 and 20 cm. The observed fluxes of 3.0, 5.3, 
8.2 and 11 mJy nearly exactly follow a power law. (By fitting 
such a law through the four data points, we find a spectral in- 
dex of -0.7 ± 0.1). A consequence of this is that there is no 
significant upper bound on the shock position Rs . (Only when 
the larger-wavelength fluxes fall below the power-law extrapo- 
lation of the short-wavelength fluxes can a constraint on Rs be 
derived - see Van Loo (2005 1.) The fits show that a large range 
of shock positions and compression ratios fit the observations. 
As for the AC 116 data, we find that the velocity jump is too 
high (Am > 500 km s ) to be compatible with hydrodynamical 
models. 

We also fitted the observations of the AB1065 dataset, 
which has detections at 3.6 and 6 cm, and a significant upper 
limit at 20 cm. Results from the fits show that the position of 
the shock is between 50 and 200 R*. The compression ratio is 
not as well constrained (a consequence of having only an upper 
limit at 20 cm): we find^ > 2. Solutions are found for all mag- 
netic field strengths B* > 10 G. Most important, however, is 
that the Au values are about the same as for the AC1 16 dataset. 

The observation AB1048 is representative of the low state. 
However, no significant constraints can derived from this ob- 
servation. The single shock velocity jump should be at least 
190 kms -1 . This value is lower than for the others because the 
general flux level for this observation is low. 

The location of the AB1065 shock is more inward than the 
AC1 16 shock. This suggests a direct explanation for the differ- 
ence in flux level between the two observations. The AB1065 
shock is closer to the star, so the intrinsic synchrotron flux has 
to pass through more free-free absorbing material, resulting in 
a lower radio flux. If this is the correct explanation for the vari- 
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ability in the observed fluxes, we should expect the variability 
to happen on the hydrodynamical time-scale, which is of the 
order of 4 days (the time needed for a shock to travel 100 R*). 
However, the AB67 1 dataset (Sect. 12. 2\ does not show any vari- 
ability. To explain why the flux remains constant for nearly one 
month, we would have to invoke the fortuitous circumstance 
of a train of shocks, all of about equal strength, occurring dur- 
ing one month. Ascribing the variability in the whole dataset 
to differences in shock positions therefore does not appear con- 
vincing. 

As an alternative cause of variability, we could consider 
variations in the magnetic field. Very little is known about mag- 
netic fields in O stars, so it is difficult to make predictions on 
what time-scale the variability would occur. One possibility is 
the rotational time-scale (which is about 1-2 weeks for these 
stars), but this is again contradicted by the AB671 observations 
(Sect. El. 

In summary, the single-star scenario can explain the spec- 
tral shape, but has considerable difficulties explaining the flux 
levels, as a high velocity jump or a very large number of 
shocks would be required. We cannot attribute the variability 
to changes in the shock structure. In the single-star scenario, 
we therefore have to surmise that the variability is caused by 
an inherent variability in the magnetic field. The period of the 
variability is discussed in Sect. 13.21 

3.2. Binary scenario 

In this section we explore the constraints the data provide as- 
suming that HD 168112 is a colliding-wind binary. The best 
studied example of such a colliding-wind binary is WR 140 
(White & Becker 1995 1 and we will use it as a guide in in- 
terpreting the present observations. The 2 cm light curve of 
WR 140 shows, during one period, a rising branch followed 
by a descending branch. At larger wavelengths, the same fea- 
ture is seen, but the peak formed by the two branches becomes 
sharper. The light curve shows good repeatability from one pe- 
riod to another. 

The intrinsic synchrotron radiation will depend on the sepa- 
ration of the two components. In an eccentric binary more syn- 
chrotron emission will be generated at periastron than at apas- 
tron. However, free-free absorption will remove some, or all, 
of the flux, depending on the position in the wind (i.e. on the 
optical depth along the ray to the observer). The observed radio 
lightcurve will therefore have a maximum at a phase different 
from periastron. Exactly when the maximum occurs depends 
on the geometry of the orbit. 

For close binaries, the colliding-wind region will be com- 
pletely inside the free-free radius, and we therefore expect to 
see no variation on the time-scale of the orbital period. This 
is indeed the case for Cyg OB 2 No. 5 (Miralles et al. 119941 . 
For HD 1681 12, however, we may assume that it is not a close 
binary, otherwise De Becker et al. (2004 1 would have found 
radial velocity variations in their optical spectra, even for an 
unfavourable orbital inclination (see also below). 

In the binary hypothesis, we therefore attribute the large 
radio flux changes in Fig. [l\ to the orbital motion. The peaks 
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Fig. 4. The 3.6 and 6 cm radio fluxes, folded with a 1.4 yr pe- 



riod. The normalised X-ray fluxes from De Becker et al. ( 2004 1 
are also plotted (bottom panel). Phase 0.0 was arbitrarily set at 
JD = 2444000.0. 



in flux occur at 1985.0, 1993.1 and 2001.8, which are ~ 8.5 yr 
apart. The true period is much shorter than this, as is shown by 
the rising and descending branches of the lightcurve. In 1984-5 
and 1992-3 we have a rising branch (seen in 2, 3.6 and 6 cm) 
lasting ~ 1 yr. A descending branch (seen in 3.6, 6 and 20 cm) 
is seen around 2002 and lasts ~ 0.45 yr. We therefore expect 
the true period to be between 1 and 2 years. 

To better determine the period we tried various periods in 
the range 1-2 yr and plotted the 3.6 and 6 cm fluxes as a func- 
tion of phase. (The 2 and 20 cm fluxes do not provide addi- 
tional information as they are strongly correlated with the 3.6 
and 6 cm fluxes.) We do not have a mathematical description 
for the expected behaviour of the fluxes as a function of time, so 
the quality of the phase locking was judged by eye. We checked 
that the lightcurve shows only one rising and one descending 
branch over the whole period and that fluxes at nearby phases 
have similar values. Trying a large number of periods we find 
that P = 1 .07, 1.12, 1 .22, 1 .4 and 1 .66 yr give good phase lock- 
ing. The result for P = 1.4 yr is shown in Fig.|4] 

Even the longest period proposed would result in an orbit 
of the companion that would remain within the free-free radius, 
at least if we assume the orbit to be circular. The orbit with a 
1 .66 yr period would have a semi-major axis of ~ 90 R* (based 
on the mass listed in Table ^for the primary and an assumed 
mass of 30 Mq for the secondary). At a radius of 90 R«, the 
radial optical depth for free-free absorption at 20 cm is ~ 2.5 



(which can be derived from Wright & Barlow ( 1975 1, using the 
data from TableQ. This would absorb a large part of the intrin- 
sic radio flux. Shorter periods would of course suffer even more 
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free-free absorption. We therefore need to assume an eccentric 
orbit, so that the colliding-wind region becomes visible near 
apastron. Note that De Becker et al. (2004) showed that the 
X-ray variability of HD 168112 can be explained by a binary 
provided the eccentricity is > 0.5. An eccentric binary with a 
1-2 yr period is also compatible with the AB67 1 fluxes being 
constant over a one month time-scale. 

X-ray data can provide some further constraints on the pe- 
riod. Similar to the radio, more X-ray flux is generated at pe- 
riastron than at apastron. For such a wide binary seen under 
low inclination, the absorption in X-rays is, however, less than 
that at radio wavelengths. As the separation between the two 
components becomes smaller, the observed X-ray flux there- 
fore increases, while the radio flux decreases (due to the in- 
creasing free-free absorption). For an eccentric binary, the pas- 
sage through periastron will be fast, so we expect a sharp peak 
in the X-ray emission around that time. 

For HD 1681 12, De Becker et al. d2004l present two XMM- 
Newton observations taken quasi-simultaneously with VLA ra- 
dio data (programmes AB 1048 and AB1065). They alsoreanal- 
ysed three older X-ray observations. If we fold all the X-ray 
fluxes with the radio periods found above, we find that only 
the 1 .4 yr period shows a reasonable behaviour for the X-ray 
fluxes. Figure|4]shows the X-ray data folded with the 1 .4 yr pe- 
riod. The steeply descending fluxes around phase 0.2-0.4 could 
be due to the fast change in separation as the secondary moves 
away from periastron in the eccentric orbit. This would suggest 
that periastron occurs at phase ~ 0.2. 

Clumping of the stellar wind material will not allow a 
shorter period, or a less eccentric orbit. This is because clump- 
ing has no effect on the optical depth, provided that we require 
the same (thermal) flux to be emitted. This can be seen from 
Abbott et al. (1981). Assuming that the density ratio of the 
low-density versus high-density material is zero, their Eq. (10) 
shows that the flux is proportional to (M 2 //) 2 ^ 3 , where / is the 
volume filling factor of the high-density material. The optical 
depth (their Eq. (9)) is proportional to M 2 / '/. If we compare a 
smooth wind with a clumped wind and require both to have the 
same flux, they need to have the same M 2 If value, which also 
means that the optical depth is the same (for details, see Van 
Loo 120051 . Clumping can also lead to the material becoming 
porous to radiation (Owocki et al. l2004> . If we would include 
the effect of porosity as well, the 20 cm optical depth at 90 R» 
would be smaller, allowing a shorter period, or a less eccentric 
orbit for the binary. 

One can speculate on the nature of the companion. If it is 
a normal star of the same luminosity class as the 05.5 III(f + ) 
primary, the spectral type would have to be later than 05.5. It 
would probably not be as late as 09 since otherwise we should 
see low excitation lines in the spectrum that would be untypical 
of the primary. If the difference in magnitude becomes larger 
than 2.5, then the companion could be of any type (a B-star, a 
pre-main-sequence star of a few solar masses, ...). A compact 
companion would produce much brighter X-ray emission than 
what is observed and the spectrum would follow a blackbody 
rather than an optically thin spectrum. Bieging et al. ( 1989) al- 
ready dismissed the compact companion hypothesis as improb- 
able on evolutionary grounds. 
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Fig. 5. The semi-amplitude K of the velocity variations as a 
function of period P, for eccentricity e - 0.5, mass ratio q = 
0.4 and inclination i = 15° (solid line), 30° (dotted), 60° (short- 
dashed) and 90° (long-dashed). The three horizontal lines at the 
10 km s level show the time ranges that have been covered by 
the available optical spectroscopy data (De Becker et al. 2004 1. 



Is the presence of a companion compatible with the fact 
that such a companion was not detected in the available opti- 
cal spectroscopy? We checked this by exploring the parameter 
space of possible periods P, mass ratios q = M<ilM\, eccen- 
tricities e and inclination angles i, to see which combinations 
give velocity changes that would remain undetected. We con- 
sider a grid of q = 0.1, 0.4, 0.7, 1.0, e = 0.0, 0.5 and i = 15°, 
30°, 60°, 90°. We also took into account the time-scales that are 
covered by the available optical data. Adopting an upper limit 
of 10 km s _1 (2 sigma) on the semi-amplitude (K) of the veloc- 
ity variations, we find that for q > 0.4, the inclination should 
be i < 15° for a period P < 2 yr. For q = 0.1, higher inclination 
angles are possible and even i - 90° cannot be excluded. Of 
course, there is always the possibility that the optical spectra 
might have been taken at the wrong phase, so that we sample 
only a small part of the total velocity variation. 

Thus, the observations are in agreement with HD 168112 
being a (low inclination) binary, with a period of less than 2 yr. 
It seems that the most likely period is ~ 1 .4 yr. The binary sce- 
nario explains the observed variability very well, as shown in 
Fig.0] While we did not model the spectral shape or flux lev- 
els, we can derive some information from the single-star syn- 
chrotron model (Sect. l3~Tl . Even though that model is not ap- 
propriate for a binary, the location and the velocity jump of the 
single shock in that model are in qualitative agreement with 
what can be expected for a colliding-wind region (e.g. WR 147, 
Dougherty et al. l2003l . 
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4. Conclusions 

In this paper we studied the available VLA and ATCA radio 
data on HD 1681 12. The archive data confirm the non-thermal 
nature and high variability of the radio emission. Variability is 
present at a time-scale of 0.45 yr or more, and the fluxes can 
differ by a factor of ~ 20. No variability was found on time- 
scales less than 0.1 yr. 

We interpreted the observations in terms of both a single- 
star and binary scenario. Assuming HD 168112 to be single, 
we applied the Van Loo et al. (20051 model to the observa- 
tions. The fact that the strong shocks (which dominate the syn- 
chrotron emission) are at different positions for different times 
of observation was suggested as an explanation for the variabil- 
ity. It fails, however, to explain why the flux remained constant 
for about one month. An even more important problem with 
the single-star scenario is that the shocks need to have veloc- 
ity jumps very much higher than those predicted by hydrody- 
namical models. As an alternative, when we forced the shocks 
to have an acceptable velocity jump, we found that too many 
shocks are required to explain the flux levels. The failure of the 
single-star scenario to give results compatible with hydrody- 
namical simulations is a major conclusion of this work. 

The binary hypothesis is very successful in explaining all 
the observational material. The non-thermal emission is due to 
the colliding winds. The variability can be explained by the 
colliding-wind region in an eccentric binary moving in and 
out of the region where the free-free absorption is important. 
The high value found for the shock velocity jump is to be ex- 
pected for a colliding-wind region. Also, the X-ray data are 
in agreement with this scenario. The non-detection by optical 
spectroscopy suggests an orbital inclination i < 15°. While the 
amount of radio data available is not large, they clearly show 
that the period should be between one and two years. In com- 
bination with the X-ray observations, the present data suggest 
that the binary period is most likely ~ 1 .4 yr. The constraint on 
the period is another major conclusion of this work, and it will 
allow a much better planning of future observing campaigns. 
While the current evidence favours the binary scenario, confir- 
mation of this and the determination of the period still need to 
be sought from future observations. 
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Appendix A: Data reduction 

A.1. VLA data 

The data reduction was done using the Astronomical Image 
Processing System (AIPS), developed by the NRAO. Technical 
details of the reduction are listed in Table lA. fl for the 2, 3.6 and 
6 cm observations and Table lA. 2l for the 20 and 90 cm observa- 
tions. Except where otherwise indicated, all observations were 
made in 2 sidebands (denoted IF1 and IF2), each of which has 
a bandwidth of 50 MHz. For the flux calibration, we followed 
the procedure outlined by Perley & Taylor (2003 ). The fluxes 
assigned to the flux calibrators are given in Table lA. 31 We then 
determined the instrumental gains from the phase calibrator ob- 
servations, interpolated them and applied them to the observa- 
tion of the source. Where needed, bad data points were flagged. 

The calibrated visibilities were assigned weights using ro- 
bust uniform weighting (Briggs 1995) and were then inverted 
by a Fourier transform to produce an image. For all reduc- 
tions, the image covers most or all of the full primary beam. 
The image was then deconvolved using the clean algorithm 
to remove the effect of the beam (point spread function). We 
stopped cleaning when the algorithm started finding about the 
same number of negative as positive components. On the obser- 
vations made with low spatial resolution (VLA configurations 
C and D), some Galactic background structure can be seen. We 
excluded the minimum amount of observational data taken on 
the shortest baselines to ensure the removal of this background. 
In Fig. IA.1I we show some typical examples of the resulting 
maps. 

We measured the fluxes on the cleaned image using the 
AIPS task jmfit to fit elliptical Gaussians to the source. The 
elliptical Gaussian was forced to have the major axis, minor 
axis and position angle equal to the beam values. For the 20 cm 
AO 16 (1984-11-27) observation, HD 168112 is not close to 
the centre of the map. In that case, the flux measurement be- 
comes more complicated. First, the sensitivity of the image de- 
creases as one moves away from the centre ("primary beam at- 
tenuation"). The images can easily be corrected for this effect 
(using the task pbcor). A second effect is beamwidth smearing 
(Bridle & Schwab 1999): an off-centre source is smeared out 
in the direction that is radial toward the centre of the image. 
Fortunately, this effect conserves the total flux. For this obser- 
vation, we measured the flux by using an elliptical Gaussian 
with all parameters left free. As a check, we used the interac- 
tive task tvstat to integrate the flux over the area occupied by 
the beamwidth-smeared source. We found that both measure- 
ments were in acceptable agreement. 

Usually, the measured fluxes are assigned an error bar cor- 
responding to the root-mean-square (RMS) noise in the map. 
The RMS error bar only covers the random sources of error. To 
estimate some of the systematic errors, we redid the reduction 
a number of times. We excluded one antenna from the data, ex- 
cluded one run (if there were at least three runs), tried natural or 
uniform weighting instead of robust uniform weighting or var- 
ied the extent of the taper. The error bar we assign to the mea- 
sured fluxes covers the range in fluxes we thus found. In many 
cases this is larger than the RMS error bar. We also measured 
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Table A.l. Reduction of the data at 2, 3.6 and 6 cm. Programme C978 is an ATCA observation, all others are VLA observations. 
Column (1) gives the programme name, (2) the date of the observation, (3) the phase calibrator (J2000 coordinates), (4) the 
distance of HD 1681 12 to the phase calibrator (degrees), (5) the integration time (in minutes) on the source, (6) the number of 
antennas that gave a usable signal, (7) the configuration the VLA or ATCA was in at the time of the observation, (8) the major 
and minor axis (both in arcsec) of the beam, (9) the position angle (degrees) of the beam, (10) the RMS in the centre of the image 
and (1 1) the measured flux (in mJy). 



(1) 


(2) 


(3) (4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


progr. 


date 


phase calibrator 


intgr. 


no. 


config 


beamsize 


PA 


RMS 


flux 






name dist. 


time 


ants. 




(arcsec 2 ) 


(deg) 


(mJy) 


(mJy) 



2 cm 






















AA29 


1984-04-04 


1733-130 


11.2 


22 


24 


C 


2.3 xl.4 


13.9 


0.20 


1.2 + 0.3 


AO 16 


1984-12-21 


1733-130 


11.2 


45 


27 


A 


0.16x0.12 


2.2 


0.13 


2.4 ± 0.2 


3.6 cm 






















BP1 


1992-05-30 


1733-130 


11.2 


43 


24 


C 


9.2 x4.8 


-88.8 


0.05 


0.30 ± 0.05 


AB671 


1993-01-21 


1811-209 


9.0 


19 


22 


A 


0.33x0.21 


25.8 


0.05 


3.42 ± 0.09 


AB671 


1993-01-24 


1811-209 


9.0 


19 


27 


A 


0.41x0.23 


38.2 


0.05 


3.51 ±0.09 


AB671 


1993-01-29 


1811-209 


9.0 


19 


26 


BnA 


0.91x0.50 


51.9 


0.05 


3.32 ± 0.09 


AB671 


1993-02-01 


1811-209 


9.0 


6 


27 


BnA 


0.82x0.40 


84.4 


0.09 


3.1 ±0.2 


AB671 


1993-02-06 


1811-209 


9.0 


19 


27 


BnA 


0.83x0.42 


68.3 


0.04 


3.30 ±0.11 


AB671 


1993-02-14 


1811-209 


9.0 


9 


27 


BnA 


0.95x0.49 


56.3 


0.09 


3.4 ± 0.2 


AB671 


combined 












0.60x0.29 


64.6 


0.02 


3.39 ± 0.07 


C978 


2001-10-11 


1832-105 


3.7 


161 


6 


EW352 


4.51x0.65 


-1.3 


0.25 


3.0 ± 0.5 


AB1048 


2002-03-24 


1832-105 


3.7 


11 


26 


A 


0.31x0.25 


11.3 


0.06 


0.36 ± 0.07 


AB1065 


2002-09-11 


1832-105 


3.7 


9 


26 


CnB 


3.0 x2.1 


47.4 


0.11 


1.85 ±0.12 


6 cm 

AA28 


1984-03-09 


1733-130 


11.2 


28 


26 


CnB 


3.7 x2.1 


-87.8 


0.07 


1.63 ±0.09 


AA29 


1984-04-04 


1733-130 


11.2 


14 


27 


C 


5.9 x3.7 


10.6 


0.11 


1.81 ±0.12 


AO 16 


1984-12-21 


1733-130 


11.2 


50 


27 


A 


0.55x0.36 


-20.1 


0.06 


5.64 ±0.13 


AB671 


1993-01-21 


1811-209 


9.0 


19 


24 


A 


0.73x0.68 


-27.2 


0.07 


4.24 ±0.11 


AB671 


1993-01-24 


1811-209 


9.0 


19 


27 


A 


0.55x0.44 


9.6 


0.06 


4.52 ±0.11 


AB671 


1993-01-29 


1811-209 


9.0 


19 


26 


BnA 


1.7 xl.l 


34.7 


0.06 


4.30 ±0.11 


AB671 


1993-02-01 


1811-209 


9.0 


6 


27 


BnA 


1.5 xO.81 


82.5 


0.09 


4.09 ±0.12 


AB671 


1993-02-06 


1811-209 


9.0 


18 


27 


BnA 


1.5 xl.O 


64.1 


0.06 


4.36 ±0.15 


AB671 


1993-02-14 


1811-209 


9.0 


10 


27 


BnA 


1.7 xO.92 


51.3 


0.08 


4.36 ±0.14 


AB671 


combined 












0.91x0.47 


59.1 


0.03 


4.34 ±0.10 


C978 


2001-10-11 


1832-105 


3.7 


161 


6 


EW352 


7.41x1.09 


-1.4 


0.19 


5.3 ± 0.45 


AB1048 


2002-03-24 


1832-105 


3.7 


11 


26 


A 


0.53x0.36 


12.3 


0.07 


0.27 ± 0.07 


AB1065 


2002-09-11 


1832-105 


3.7 


8 


26 


CnB 


4.8 x2.0 


57.4 


0.13 


2.32 ±0.15 



BPl 1992-05-30 3.6 cm AB671 combined 3.6 cm AB671 combined 6 cm 




18 h 18 m 42 s 40 s 18 h 18 m 41.0 s 40.8 s 18 h 18 m 41.0 s 40.8 s 

RA (J2000) RA (J2000) RA (JS000) 

-0.1,0.1,0.2,0.3 mJy -0.05,0.05,0.1,0.2,0.5,1.0,2.0,3.0 mJy -0.05,0.05,0.1,0.2,0.5,1.0,2.0,3.0,4.0 mJy 

Fig. A.l. Some examples of maps of VLA archive data at 3.6 and 6 cm. The values of the contour levels are listed at the bottom 
of each panel. The negative contour is given by the dashed line. The first positive contour is at about twice the RMS noise in the 
map. The beam is shown in the upper left corner. The cross indicates the optical ICRS 2000.0 position (from SIMBAD). Because 
self-calibration loses absolute positional information, this is not shown in the middle and right panels. Note the different spatial 
scales and different contour levels of the maps. 
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Table A.2. Reduction of the data at 13, 20 and 90 cm. Programmes C599 and C978 are ATCA observations, all others are VLA 
observations. Columns are as in Table lA.fl note that column (10) gives the RMS noise for the standard reduction in the centre 
of the image, while upper limits in column (11) are based on the highest RMS found in the search for systematic effects (see 
Appendix lA.H and are corrected for the primary beam effect and bandwidth smearing, if not in the centre of the image. Column 
(12) refers to the notes and, for the VLA observations, lists the width of the sidebands if different from 2 x 50 MHz. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


nropr 


date 


phase calibrator 


int?r 


no. 


config 


beamsize 


PA 


RMS 


flux 


notes 






name 


dist. 


time 


ants. 




(arcsec") 


(dea) 


(mlv) 






13 cm 
























C978 


2001-10-1 1 


1832-105 


3.7 


120 


6 


EW352 


15. x 1.9 


-1.3 


0.26 


8.2 ± 0.6 


X 


20 cm 
























BOWE 


1982-01-31 


1743-038 


11.9 


13 




A 


4 2x 1.2 


75.2 


7.1 


< 100. 


R 1 x 1 562^ 


BAUD 


1982-02-28 


1743-038 


11.9 


33 


13 


A 


1.4x 1.0 


5.5 


11. 


< 55. 


B 1 x 1 5625 


AO 16 


1QR4-1 1-27 

X J U i XX — 1 


1733-130 


1 1.2 


28 


25 


A 


1.7x 1.2 


9.4 


3.6 


11+7 


B 


AO 16 


1QR4- 12-21 

X J U i X — 1 


1733-130 


1 1.2 


10 


27 


A 


1.7x 1.2 


3.4 


0.14 


8.9 ± 0.2 




AG 163 


1984-12-24 


1743-038 


1 1.9 


1 


27 


A 


1.8x 1.2 


-9.8 


0.47 


< 90. 


B,2 x 25 


AO 16 


1985-02-16 


1733-130 


1 1.2 


15 


25 


A 


1.6x 1.2 


-2.4 


0.12 


< 12. 


B 


AB531 


1 989-05-01 


1834-126 


3.9 


2 


26 


B 


6.3x 4.1 


17.9 


1.0 


< 26. 


B,l x 50 


VP51 


1989-06-19 


1733-130 


1 1.2 


29 


25 


c 


15. xlO. 


-14.4 


0.31 


< 1.8 


X 


AG290 


1989-07-20 


ion Tin 

1833-210 


9.6 


o 


27 


c 


22. xll. 


-21.1 


1.4 


< 60. 


B 


VP91 


1989-1 1-10 

i yuy i± xv/ 


1733-130 


1 1.2 


18 


27 


Q 


43. x27. 


-12.2 


1.7 


< 6. 


x 


BP1 


1992-05-30 


1011 901 

1711 1 


1 4 


40 


25 


c 


20. xl5. 


-70.5 


0.80 


< 3.8 




TSTOB 


1993-07-05 


1822-096 


2.6 


14 


27 


c 


67. xl2. 


50.2 


0.84 


< 3.7 




AS534 


1994-05-26 


1822-096 


2.6 


5 


24 


BnA 


4.1x 2.3 


-79.5 


1.2 


< 8. 


B,X,1 x 3.125 


AS534 


1994-06-01 


1822-096 


2.6 


4 


25 


BnA 


4.9x 2.1 


56.8 


1.6 


< 90. 


B,X,1 x 3.125 


AC308 


1996-06-09 


1833-210 


9.6 


1 


27 


DnC 


38. x21. 


81.1 


2.4 


< 19. 


B 


AC308 


1996-06-16 


1833-210 


9.6 


1 


27 


DnC 


40. x20. 


71.7 


1.6 


< 18. 


B 


C599 


1997-02-25 


1934-638 


54.0 


89 


6 


6A 


32. x 4.4 


3.2 


0.52 


<23. 


B 


AC308 


1997-10-16 


1833-210 


9.6 


1 


27 


DnC 


38. x20. 


84.0 


2.2 


< 16. 


B 


C978 


2001-10-11 


1832-105 


3.7 


120 


6 


EW352 


27. x 3.4 


1.2 


0.42 


11 ± 2 


X 


AB1048 


2002-03-24 


1834-126 


3.9 


5 


26 


A 


2.0x 1.3 


19.0 


0.17 


<0.7 




AB1065 


2002-09-11 


1834-126 


3.9 


5 


26 


CnB 


7.1x 3.3 


40.0 


0.68 


< 2.8 


I,A = 18 cm 


AB1065 


2002-09-11 


1834-126 


3.9 


5 


21 


CnB 


8.1x 3.8 


38.9 


0.37 


< 1.5 


A = 20 cm 


90 cm 
























AH299 


1988-06-21 


1829+487 


60.9 


47 


21 


DnC 


125. x72. 


-81.0 


57. 


< 340. 


B,2x 3.125 


AH299 


1989-05-28 


1829+487 


60.9 


57 


25 


CnB 


50. x37. 


-65.7 


13. 


<90. 


B,2x 3.125 



Notes: 

B Observation is not centered on HD 1681 12. The flux (or upper limit) has been corrected for that. 

No correction for beamwidth smearing needs to be applied to C599 (see Sect. IA.3l . 
X The phase calibrator is of low quality (Perley & Taylor[3UU3J for VP5 1 , VP9 1 , AS534 and C978. 
I Interference is high in the AB1065 18 cm image. 



the flux by determining the maximum intensity of the point 
source, instead of fitting an elliptical Gaussian. In most cases, 
both measurements agree within the error bars. The worst non- 
agreement is the 6 cm AC1 16 observation, where the maximum 
intensity is -2.5 <x below the integrated value. Lesser effects are 
shown by AB1065 - 3.6 cm (+1.2 cr), AA28 - 6 cm (-1.4 cr) 
andAB671 1993-01-24- 6 cm (-1.5 cr). 

The final error bar was determined by adding (in a RMS 
sense) a 2 % calibration error for the 3.6, 6, 20 and 90 cm ob- 
servations and a 5 % error for the 2 cm. These represent the 
uncertainty in the absolute calibration (Perley & Tavlor l2T)03> . 

For those observations where HD 168 1 12 was not detected, 
we assigned an upper limit of 3 times the RMS noise mea- 
sured in the map (in a small box around the source). To take 
into account some of the systematic effects, we also performed 



the various reductions described above and used the highest 
RMS thus found to define the upper limit. When HD 168112 
is offset from the centre, the RMS must be measured on the 
image corrected for primary beam attenuation. To correct for 
the beamwidth smearing, we had to multiply the RMS by the 
number of beams that a detectable beamwidth-smeared source 
at this position would cover. From Bridle & Schwab (1999 1 
we find that the number of beams can be approximated by 
+/3 2 , where: 

bandwidth (MHz) offset (arcsec) 

/3 = x — . (A.l) 

observing frequency (MHz) beamwidth (arcsec) 

In some cases HD 168112 is within the primary beam on 
two images centred on different targets. In that case we only 
list the one closest to HD 1681 12 because this gives the small- 
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Table A.3. Fluxes of the flux calibrators used in the re- 
duction of the archive data. Programmes C599 and C978 
are ATCA observations, all others are VLA observations. 
The VLA flux scale is based on the 1995.2 coefficients 
(Perley & Taylor |2"UU51 . The ATCA flux scale is given by 
Reynolds (1993). 



calib. 


freq. 


flux 


programmes 




(GHz) 


(Jy) 




2 cm 








3C286 


14.915 


3.460 


AA29,AC116 




14.965 


3.452 


AA29,AC116 


3.6 cm 








3C286 


8.415 


5.224 


AB671,BP1 




8.465 


5.203 


AB671 


3C48 


8.435 


3.164 


AB 1048 AB 1065 




8.485 


3.145 


AB 1048 AB 1065 


1934-638 


8.640 


2.841 


C978 


6 cm 








3C286 


4.835 


7.510 


AA28,AA29, AC 1 16 AB67 1 




4.885 


7.462 


AA28,AA29,AC1 16,AB671 


3C48 


4.835 


5.459 


AB 1048 AB 1065 




4.885 


5.405 


AB 1048 AB 1065 


1934-638 


4.800 


5.828 


C978 


13 cm 








1934-638 


2.496 


11.14 


C978 


20 cm 








3C286 


1.365 


15.01 


AC308 




1.435 


14.65 


AC308 




1.452 


14.57 


AG 163 




1.465 


14.51 


AC116,AG290,AB531 




1.502 


14.34 


AG 163 




1.515 


14.28 


AC116AG290 




1.612 


13.85 


BAUD,BOWEAS534 




1.665 


13.63 


VP91,BP1 


3C48 


1.385 


16.20 


AB 1048, AB 1065,TSTOB 




1.465 


15.49 


AB1048.AB1065 




1.635 


14.19 


AB1065 




1.665 


13.98 


AB1065 


1733-130 


1.665 


5.200 


VP51 a 


1934-638 


1.384 


14.94 


C599,C978 


90 cm 








3C286 


0.328 


25.90 


AH299 




0.333 


25.77 


AH299 



a for VP51, no flux calibrator was available, so we calibrated the 
flux on the phase calibrator, to which we assigned the value listed 
in Perley & Taylor <2003l . 



est error bar. We did check that the flux of the other image 
is compatible with the flux listed. The following observations 
were not listed because they are of such low quality that they 
do not even provide a significant upper limit: BAUD (1980- 
07-12; 20 cm), FIX (1982-05-22; 20 cm), AC161 (1986-05- 
11; 90 cm), AH250 (1986-12-27 and 1987-03-19; 90 cm), 
AK162 (1987-03-20; 90 cm), AH299 (1988-03-24, 1988-05-02 
and 1988-09-23; 90 cm), AT143 (1993-06-11; 20 cm), AL341 
(1994-09-13; 90 cm) and AG626 (2002-08-14, 2002-11-23, 
2002-12-07, 2003-02-23 and 2003-02-24; 90 cm). 



A complication arose in the reduction of the TSTOB data. 
Comparing the fluxes of other sources in the field around 
HD 168112 across the various programmes showed that the 
TSTOB values are too low by a factor of ~ 2. Fortunately, 
the TSTOB programme was immediately followed by another 
programme (AW362) that observed the same flux calibrator 
at the same wavelength. Using the flux calibrator from that 
programme in the reduction, the TSTOB fluxes become high 
enough to be compatible with the others. The results in Fig.^ 
and Table lA.2l are based on this improved reduction. 

A comparison with the literature (Bieging et al. 1989 
Phillips & Titus [T990I De Becker et al. 120041 shows accept- 
able agreement in most cases. There is a tendency for our error 
bars (or upper limit in the case of Phillips & Titus) to be larger: 
this is due to our inclusion of systematic effects. Note also that 
the AB1048 and AB1065 fluxes are somewhat different from 
those published in De Becker et al. (2004). This is because, in 
the present paper, we followed the recommendation by Perley 
& Taylor ( 2003 1 to assign a small, but non-zero, weight (wtuv 
= 0.1) to the flux calibrator visibilities outside the usable range. 
There are two exceptions to the acceptable agreement. For the 
2 cm AC1 16 (1984-12-21) observation, we have 2.4 + 0.2 mJy, 
while Bieging et al. find 1.3 ±0.1 mJy. They note an apparent 
increase of flux at the largest baselines. Our reduction of the 
data does not show such a problem, giving a higher confidence 
in our result. For the 6 cm AA28 (1984-03-09) observation we 
found 1.63 + 0.09 mJy, while Bieging et al. list 1.3 + 0.1 mJy. 
We could not find an explanation for this discrepancy. 

A.2. Self-calibration 

The AB671 dataset has a number of problems. The phase cal- 
ibrator used has a low 3.6 cm flux, making it unsuitable for 
phase calibration. The time between two phase calibrations ex- 
ceeds 30 min in 5 of the 6 observations, which is too long to 
follow the phase changes due to the atmosphere, especially in 
these high-resolution observations. When we applied our stan- 
dard reduction, we found considerable variability in the fluxes. 
However, these flux variations are highly correlated with the 
flux variations of HD 167971, which was observed as part of 
the same observing programme. At least a large part of the vari- 
ability must therefore be due to atmospheric effects that were 
not correctly removed in the standard reduction. 

We therefore decided to apply phase self-calibration (see, 
e.g. Cornwell & Fomalont 1999) to this set of observations. 
This technique is suitable for the "simple" field here, consist- 
ing of two point sources on an empty background. It allows 
us to follow the atmospheric phase changes on a shorter time- 
scale than the time between two phase calibrator observations. 
From the known sensitivity of the VLA and the estimated flux 
of HD 168112, we find that a 3 min time-scale is appropri- 
ate, provided that we combine both sidebands and all Stokes 
polarizations and decrease the requested signal-to-noise ratio 
of the gain solutions to 2 (Ulvestad|2002i. Ideally, we would 
prefer a 10 sec time-scale (which is the duration of a single 
integration), so we could follow the phase changes even bet- 
ter. Unfortunately, insufficient signal would then be available 
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C978 



3.6 cm 




C978 2001-10-11 



6 cm 



-12 u 06'30" - 




C978 



4L0 b 40. 6 b 

RA (J2000) 
-0.5,0.5,1,2 mjy 

2001-10-11 



13 cm 



C978 



18"l8 jli 41.0 b 40.0 b 
RA (J2000) 
-0.5,0.5,1,2,3,4,5 mjy 

2001-10-11 20 cm 




-12"07' 0" - 



18"18 J11 42.0 S 40.0° 
RA (J2000) 
-0.5,0.5,1,2,4,6 mjy 



18"l8 111 42.0 b 40.0 s 
RA (J2000) 
-1,1,2,4,6,8,10 mjy 



Fig.A.2. Maps of ATCA archive data where HD 168112 was 
detected. The values of the contour levels are listed at the 
bottom of each panel. The negative contour is given by the 
dashed line. The first positive contour is at about twice the 
RMS noise in the map. The beam is shown in the upper left 
corner. The cross indicates the optical ICRS 2000.0 position 
(from SIMBAD). Note the different spatial scales and different 
contour levels of the maps. 



for the calibration. We applied 5 iterations of phase-only cal- 
ibration. No amplitude self-calibration was done because am- 
plitude errors are relatively unimportant on the dynamic ranges 
of the images we have. The results presented in Figs.^andEI 
and Table lA. fl are those based on the self-calibration reduction. 

To make sure that the self-calibration technique does not 
introduce systematic effects, we also applied it to those cases 
where it is not really needed. These tests show that self- 
calibration and the standard reduction give the same results 
within a few percent, provided that the self-calibration is care- 
fully applied (i.e. not on a complicated background, no inter- 
ference present in the image and a time-scale appropriate for 
the flux). A side effect of self-calibration is that the absolute 
positional information is lost. 

A3. ATCA data 

The ATCA data were reduced in Miriad following the user 
guide (Sault & Killeen fT999b . We did not reduce the 13 cm 
data of the C599 programme because HD 1681 12 is well out- 
side the primary beam. The data were read into Miriad and cor- 
rected for self-interference of the array and for the phase differ- 
ence between the X and Y channels. After reading in, the data 
consist of 13 frequency channels, each of 8 MHz width. The 
reduction then proceeds in a way similar to the VLA reduc- 
tion ( Appendix IA. 11 1. The use of multiple frequency channels 
requires us to do a bandpass calibration on the phase calibrator. 



Fluxes assigned to the flux calibrator are listed in Table lA. 31 
The programme C978 uses the phase calibrator 1832-105 
(J2000), which is not recommended at 20 cm (according to the 
VLA Calibrator Manual, Perley & TaylorFJUUS}. For the 1 3 and 
20 cm observations, we therefore made a map around this cal- 
ibrator. To remove structure in the background, we had to flag 
the shortest baselines. In our data, the phase calibrator turns out 
to be a point source (with a flux of 1 .04 Jy at 13 cm and 0.94 Jy 
at 20 cm). There are no other sources in the field that have a 
flux higher than 0.01 Jy. Thus, we consider 1832-105 to be an 
acceptable calibrator for our data. 

For the target observation, we also had to exclude data from 
the shortest baselines, to remove background structure in the 
resulting image. In the inversion and cleaning step, we applied 
multi-frequency synthesis (MFS), which compensates for the 
spectral index of the source across the bandwidth. An inter- 
esting consequence of MFS is that we do not have to correct 
for bandwidth smearing for those targets that are offset from 
the centre (Cotton 1999 1. The cleaning was applied in a box 
around the target. At larger wavelengths more sources are vis- 
ible on the image, and more boxes were therefore used. We 
stopped cleaning when the algorithm started finding about the 
same number of negative as positive components. Applying 
this criterion, we arrive close to the theoretical RMS, which 
is the noise calculated taking only the system temperature of 
the front-end receiver into account, not the calibration errors, 
side-lobes or any other instrumental effects. Only for the 20 cm 
observation do we end up with an RMS level of about twice 
the theoretical RMS. We attribute this to the more complicated 
background at this wavelength. Results are listed in Tables |A~T1 
and lA.2l These results include the flux calibration uncertainty 
of 2 % and the estimated effect of systematic errors (derived in 
a similar way as for the VLA data). The detections are plotted 
inFig.lA2l 
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